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ABSTRACT 

We determine X-ray temperatures and abundances for elliptical galaxies 
drawn from a complete, optically selected sample. The optical magnitude- 
limited sample consists of 43 galaxies, complete to a (corrected) B magnitude 
of = 11.36. Of these, 30 have enough X-ray spectral counts to allow 
temperature determinations. We find that the temperatures of the X-ray 
emitting gas in these ellipticals are correlated with the central stellar velocity 
dispersions: T oc o"^'^^; this is a shallower trend than the simple thermal relation 
T oc cr^. The diffuse gas is substantially hotter than the kinetic temperature 
of the luminous stars, since kT > ^nipa^ for all galaxies with measurable 
temperatures. This strongly indicates that dark matter halos are characterized 
by velocity dispersions which exceed those of the luminous stars. We see no 
evidence of emission from X-ray binaries becoming progressively more dominant 
in lower luminosity ellipticals. We find that so-called "supersoft" sources adhere 
to the observed kT — a relation, so they are no softer than expected for their 
velocity dispersions. We also find that ISM temperatures and abundances are 
correlated, with the gas in hotter systems being more enriched than in cooler 
galaxies. However, no correlation is found between gaseous abundances and 
stellar abundances, as inferred from Mg2 indices. 



Subject headings: galaxies: abundances — galaxies: elliptical — galaxies: ISM 
— X-ray: galaxies 
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1. Introduction 

Elliptical galaxies have been known as X-ray sources since several were detected in the 
Virgo cluster by the Einstein Observatory (Forman et al. 1979). Long and Van Speybroeck 
(1983) compiled many more Einstein detections of ellipticals and Nulsen, Stewart and 
Fabian (1984) showed that diffuse gas is their dominate source of X-ray emission. Analysis 
of Einstein Observatory IPC spectra shows that the emission from luminous ellipticals can 
be fit by thermal plasma models with fairly low temperatures, 0.5 - 2.0 keV (Forman, Jones, 
& Tucker 1985; Trinchieri, Fabbiano & Canizares 1986). More recent Ginga, ROSAT and 
ASCA observations of individual galaxies confirm these characteristic temperatures (cf. 
Awaki et al.l991; Forman et al. 1993; Fabbiano, Kim & Trinchieri 1994; Awaki et al.l994; 
Matsushita et al.l994). The emission from lower luminosity ellipticals is apparently more 
complex, with evidence for both harder and "supersoft" emission according to Fabbiano, 
Kim & Trinchieri (1992); these authors attribute the hard component to stellar X-ray 
binaries and suggest the supersoft emission may be from M stars, ultrasoft X-ray binaries 
and/or relatively cool ISM. 

To date, there has been no study of the X-ray properties of a complete sample of 
normal elliptical galaxies. The IPC detector on the Einstein Observatory was used to 
observe ~ 200 normal galaxies, including many ellipticals (see the review by Fabbiano 
1989), but these data allow only small {N ^ 14) complete samples of eUipticals to be 
constructed. Furthermore, the limited energy resolution of the IPC prevented accurate 
temperature determinations. We have therefore turned to the longer exposures of ROSAT 
pointed observations to construct the largest possible complete sample of accurate X-ray 
luminosities and temperatures for elliptical galaxies. In this letter we present new gas 
temperature and abundance determinations from ROSAT PSPC spectra of a subset of this 
complete sample and explore their relation to optical and dynamical properties of the host 
galaxies. We will present detailed analyses of their X-ray luminosities in a forthcoming 
paper (Davis & White 1996). 

2. Sample Selection 

We used the HE AS ARC BROWSE software to select our optical magnitude-limited 
sample from the Third Revised Catalog of Galaxies (de Vaucouleurs et al. 1991; hereafter 
RC3), using its corrected B magnitudes. We limited our sample to eUiptical galaxies (de 
Vaucouleurs T-types -4, -5 & -6). We then used BROWSE to cross-correlate this sample 
with the ROSAT and Einstein pointed observations. Neglecting galaxies in the zone-of 
avoidance (taken here to be |6| < 20°), we found that 38 of the 43 optically brightest 
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ellipticals have ROSAT PSPC data, and 4 of the remaining 5 have Einstein IPC data 
(the remaining galaxy, NGC 7507, will soon have ROSAT HRI data available, but this 
has no spectral information); since ROSAT PSPC data have better spectral resolution 
than Einstein IPC data, they are preferable when there is a choice. The sample of 
43 optically-selected ellipticals is complete to a limiting (corrected) blue magnitude of 

= 11.36 and includes four Local Group dwarfs. The X-ray counterpart to our optical 
sample is thus 98% complete, with X-ray data for 42 of the 43 elhpticals. If we press 
to fainter optical magnitudes, the marginal completeness in the X-ray data deteriorates 
at the 50% level for the next 10 or so optically brightest galaxies. In this letter we 
present a subset of 30 galaxies from this complete sample which have sufficient spectral 
counts to accurately determine their X-ray temperatures; 27 of these galaxies also have 
well-constrained abundances. 



3. X-Ray Data 

In order to tie the X-ray properties of ellipticals to their optical properties in a 
uniform manner, we let the optical photometry dictate the spatial size of the X-ray spectral 
extraction regions. For each galaxy in our optically selected sample, wc extract X-ray 
spectra from within 6re, where Tg is the galaxy's optical effective radius (derived from the 
RC3). Due to X-ray scattering in the ROSAT mirror assembly, the minimum suggested 
extraction radius for point sources is 1.5', which is the 99.5% encircled energy radius for a 1 
keV spectrum. Our minimum metric extraction radius must be Sr^ to meet this requirement 
with the current sample. The purpose of our larger metric extraction radius of Org is to 
provide more global X-ray properties, for better comparison with global optical quantities 
such as magnitudes (6re encloses 91% of a de Vaucouleurs r^/^ profile). 

For the PSPC data, source spectra were selected from observation times when the 
Master Veto Rate (MVR) was within recommended bounds: 20 < MVR < 170. The 
background for each spectrum was obtained from an annulus around the galaxy extending 
from 6re to 12re from its center. All detectable point sources were excluded from both 
source and background spectra; excluded points in the source regions tended to be in the 
outskirts so they are unlikely to be X-ray binaries. Source and background spectra were 
also corrected for vignetting effects and for the residual particle background (in the manner 
prescribed by Snowden et al. 1994), which is proportional to the MVR. To ensure correct 
Poisson errors were used throughout, we modified the appropriate data header as outlined 
in the PROS QPSPEC documentation. Finally, the extracted spectra were rebinned so 
that each channel has a minimum of 25 counts, and systematic errors of 1% were included. 
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4. Analysis 

We used XSPEC 9.0 software (Schafer et al. 1990) to fit various models to the extracted 
spectra. An isothermal Raymond- Smith plasma model was our benchmark model, with 
the temperature and abundance allowed to vary in the fits. We also included a variable 
absorption component due to the column density of Galactic hydrogen in the line-of- sight. 
The redshift of the model spectrum was fixed to that corresponding to the hehocentric 
velocity of the galaxy given in the RC3. The extracted spectra were fit between ~0.2 and 
2.0 keV, the exact energy boundaries being set by the channel grouping. Once a minimum 
in was found, the 90% confidence errors were determined for the free parameters. Most 
of the galaxies have sufficient counts so that the temperature, abundance and absorption 
could all be fit. For a few galaxies, however, the number of free parameters had to be 
reduced to allow reasonable errors to be determined. For such galaxies we first fixed the 
abundance to 0.2 solar; when necessary, we also fixed the absorption to the Galactic value 
(Stark et al. 1992). The Raymond-Smith thermal model provided formally excellent fits 
to all but two galaxies, with per degree of freedom u oi xl 1- The two exceptions 
(NGC 4486 & IC 1459) are known to have AGN, so we fit composite models consisting of a 
Raymond-Smith isothermal plus a power-law spectral component; these composite spectra 
achieved xl ^ ^- The detailed spectral fitting results are shown in Table 1, where the 
indicated errors are 90% confidence limits. 



4.1. kT—a Correlation 

For the 30 galaxies which have enough counts to allow accurate temperature 
determinations. Fig. 1 shows a logarithmic plot of their temperatures kT against a, the 
central line-of-sight stellar velocity dispersions within the galaxies (from Fabcr et al. 1989 
or Dressier, Fabcr & Burstein 1991). Characteristics of the best fit linear regression to 
these data arc shown in the first line of Table 2. We found that three galaxies (M32, 
NGC 1052 and NGC 2768) have temperatures of ~ 3 — 4 keV, which are much greater 
than those of other galaxies with similar velocity dispersions. These galaxies are likely to 
have their X-ray emission dominated by X-ray binaries or AGN (cf. Eskridge, White & 
Davis 1996; Ho, Fihpenko & Sargent 1995), rather than by the cooler diffuse gas expected 
from the accumulation of stellar ejecta. Neglecting these three galaxies, the best fit linear 
regression yields kT oc (7i-45='=o.2o ^^le temperature dispersion is only 0.10 dex (see 
line 2 of Table 2); the regression is shown as a solid line in Figure 1. This correlation is 
shallower than what would obtain if the gas and stars had the same kinetic temperature: 
kT = nrripa^; this latter relation is also shown in Fig. 1 as a dashed line. The bulk of our 
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temperature measurements indicate that warm diffuse gas is being detected, even down to 

low luminosities and velocity dispersions where some investigators expect such gas to be 
swept by winds (e.g. David, Forman & Jones 1991; Ciotti et al. 1991). Fig. 1 also shows 
that the global gas temperatures are hotter than the kinetic temperatures of luminous stars 
in the galaxies, since kT > iimpa"^ in all the galaxies. 

4.2. Temperature— Abundance Correlation 

For those 27 galaxies with enough counts for elemental abundances to be constrained, 
we find that the measured abundances (driven by ~1 keV Fe L emission in ROSAT 
PSPC spectra) are correlated with the gas temperatures: hotter gas tends to have higher 
abundances (see Fig. 2). Line 3 of Table 2 shows the regression results for all the data. 
When we exclude M32 and NGC 1052 from the regression (for being hot point sources, 
rather than diffuse gas), we find the abundance A oc j'2.44±o.69 ^ggg 4 of Table 2), with 
a dispersion of 0.48 dex in A. This latter correlation is shown as a dotted line in Fig. 2 
(where M32 and NGC 1052 are also indicated). We have checked our abundances against 
the several pubhshed ASCA values and find that in most cases they agree within the 90% 
confidence level and in all cases they agree within the 99% confidence level. Our PSPC 
temperatures are also consistent with the ASCA results. We find no correlation between 
our measured abundances and the stellar abundances of the host galaxies, as inferred from 
Mg2 indices (7S) (see lines 5-6 of Table 2); this is consistent with the ASCA results of 
Loewenstein et al. (1994) and Mushotzky et al. (1994). 

5. Summciry and Discussion 

In this letter we have shown that gas temperatures are correlated with stellar velocity 
dispersions, T oc o"^'^^, and are everywhere hotter than the stellar kinetic temperatures 
[kT > umpa^). There are several possible mechanisms for the gas acquiring more specific 
energy than the stars. One possibility is that supernovae have heated the gas, but this 
is inconsistent with the generally low metal abundances being found in such galaxies 
(Matsushita et al. 1994; Loewenstein et al. 1994; and §4.2 of this paper). A more 
likely possibility is that the luminous parts of the galaxies are embedded in dark matter 
halos which are dynamically hotter than the stars (Fabian et al. 1986). As stellar ejecta 
accumulates and flows inward in cooling flows, movement through the gravitational 
potential heats the gas. If dark matter halos are characterized by velocity dispersions that 
are ^ 50% greater than those of the luminous stars, then gravitational heating can raise the 



-6- 



gas temperatures to the observed values. 

In the largest study of X-ray spectral properties of galaxies to date, Kim, Fabbiano & 
Trinchieri (1992, hereafter KFT) found that the X-ray emission from eUiptical galaxies is 
more comphcated than the simple thermal plasma found by Forman et al. (1985). KFT 
found for galaxies with high X-ray to optical flux ratios (Lx/Lb) that their spectra are 
indeed well described by Raymond-Smith thermal plasma models with temperatures of 
~ 1 keV and solar abundances. But for galaxies with intermediate values of Lx/Lb, KFT 
found spectral hardening. At even lower values of Lx/Lb, the spectra require a a very cool 
"supcrsoft" 0.2 kcV component in addition to the hard (~5 keV) component. Although a 
few galaxies in our sample have unusually hard emission, we find that the correlation of 
temperature with velocity dispersion persists down to kT 0.2 keV, Lb ~ 4.7 x lO^h'^ 
Lq, (T 120 km s"^. Our sample includes some of the galaxies deemed to have "supersoft" 
spectra by KFT, but we find that these sources are simply on the cool end of a continuous 
distribution, with temperatures consistent with expectations from ISM in galaxies with low 
velocity dispersions (hence, shallower gravitational potentials). 

While a few of the galaxies in our sample exhibit hard emission (due to X-ray binaries, 
AGN, or both) we find no general hardening of the X-ray emission as we go to lower 
luminosities and velocity dispersions, as might be expected if low-luminosity ellipticals 
are blowing winds. KFT find harder spectra for a subset of their sample (their "Group 
2"), which they attribute to X-ray binaries, but our PSPC data for 4 of 8 galaxies in this 
subsample (after removing NGC 1052, a Seyfert 2, which may have biased the results of 
KFT) show no evidence for having harder spectra attributable to X-ray binaries. Thus, 
we detect gaseous X-ray emission without contamination by X-ray binaries down to 
Lb ~ 4.7 X lO^h"'^ Lq, Lx ~ lO^^h"^ erg s~^ and a 120 km s~^. This is contrary to the 
predictions of some theoretical models for the evolution of hot gas in ellipticals, in which 
the threshold for the onset of supernovae-driven winds occurs at significantly higher optical 
luminosities (Giotti et al. 1991; David et al. 1991). Such winds drive out the gas, leaving 
only X-ray binaries as significant sources of X-ray emission. M32 may be an example of 
such a wind-blowing galaxy, since it exhibits hard X-ray emission and is ten times less 
luminous (both optically and in X-rays) than the next less luminous galaxy in our sample. 

Finally, we find that abundances are correlated with global temperatures: A oc T"^'^. 
This may refiect an underlying trend in observed stellar abundances (more luminous galaxies 
are more metal rich), but we find no detailed correlation between our gaseous abundances 
and the stellar metallicity index Mg2. The Mg2 index is measured within the inner few kpc 
of galaxies, so it may not be representative of the galaxy out to our extraction radius of 
6re. For example, Davies, Sadler & Peletier (1993) find that the Mg2 index drops by ~25% 
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from the center to a radius of ~ Ir^. Analysis of ASCA observations of the eUiptical NGC 
4636 shows there is a gaseous abundance gradient evident from spatially resolved X-ray 
spectroscopy, a gradient which persists to at least Srg (Mushotzky et al. 1994). Our global, 
emission-weighted measures of abundances obscure such structure. The lack of correlation 
between the stellar and gaseous abundances may alternatively indicate that more luminous 
galaxies have larger supernova rates than less luminous galaxies, but not so large as to 
unbind the gas. 

In a subsequent, more detailed analysis of this complete sample of ellipticals, we will 
present their X-ray luminosities, their luminosity function, discuss the Local Group dwarf 
galaxies, and assess additional correlations and implications. 

This research made use of the HEASARC, NED, and SkyView databases. We would 
like to thank Keith Arnaud for his help with the IPC data. This work was partially 
supported by the NSF and the State of Alabama through EPSCoR grant EHR-9108761. 
REW also acknowledges partial support from NASA grants NAG 5-1718 and NAG 5-1973. 
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Table 2. Correlations 
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^Probability that there is no correlation, using Kendall's r test; dexponents are given as 
exponents. 



Figure Captions 

Fig. 1. — The kT — a correlation. The solid line is the fit to the data excluding the three 
galaxies with kT > 2 keV (M32, NGC 1052, and NGC 2768). The dashed line shows the 
relation expected if the gas temperature and stellar kinetic temperatures were the same. The 
datum with the largest error bars (at log a ^ 2.42) is NGC 4105, the only galaxy in the plot 
where the temperature was determined using IPC data. 

Fig. 2. — The abundance— temperature correlation. The dashed line is the fit to the sample 
excluding the high temperature galaxies M32 and NGC 1052 (seen at right). Kendall's r 
indicates the correlation is significant at the 99.84% confidence level. 
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